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1. Introduction 

In previous studies [1,2] on the vacuum dynamics of pure non-abelian gauge theories we found 
that the deconfinement temperature depends on the strength of an external abelian chromomagnetic 
field. We also verified that the same effect is not present in the case of abelian gauge theories, sug- 
gesting that it could be linked to the non-abelian nature of the gauge group. In our opinion the 
dependence of the deconfinement temperature on applied external fields is a consequence of the 
dynamics underlying color confinement and therefore, apart from possible phenomenological im- 
plications, such an effect could shed light on confinement/deconfinement mechanisms. On these 
premises it is important to test if the effect continues to hold and how it qualitatively changes when 
switching on fermionic degrees of freedom and on that ground to investigate the dependence of 
the deconfinement temperature on the strength of an external abelian chromomagnetic field in the 
case of full QCD with two flavors. Besides another relevant issue regards the relation between 
deconfinement and chiral symmetry restoration. As it is well known, the two phenomena appear to 
be coincident in ordinary QCD, while they are not so in different theories (like QCD with adjoint 
fermions [3-5]). Even if a simple explanation of this fact is still lacking it could be strictly con- 
nected to the very dynamics of color confinement. A contribution towards a clear understanding of 
this phenomenon could be to study whether it is stable against the variation of external parameters. 
In the present we investigated whether the deconfinement temperature depends on the strength of an 
external abelian chromomagnetic field and whether deconfinement and chiral symmetry restoration 
for QCD with two flavors coincide even in presence of a constant chromomagnetic field. 



2. A constant chromomagnetic field on the lattice 

QCD dynamics in presence of a time-independent background field at finite temperature can 
be studied [6-8] by means of the free energy functional 
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with 3fr the thermal partition functional 
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The functional integration is performed over the lattice links, but constraining the spatial links 
belonging to a given time slice (say x t = 0) to be 

U k (x,x t =0)=ur t (x), (£=1,2,3), (2.3) 

L^ xt (x) being the elementary parallel transports corresponding to the external continuum gauge 
potential. In the previous equations A ext (x) is the continuum gauge potential of the external static 
background field, Sw the standard pure gauge Wilson action, Sf is the fermion action and M is the 
fermionic matrix. The spatial links are constrained to the values corresponding to the lattice version 



2 



Perturbing QCD with external fields 



Leonardo Cosmai 



of the external background field, whereas the fermionic fields are not constrained. In the case of a 
static background field which does not vanish at infinity we must also impose that, for each time 
slice x t ^ 0, spatial links exiting from sites belonging to the spatial boundaries are fixed according 
to eq. (23). In the continuum this last condition amounts to the requirement that fluctuations over 
the background field vanish at infinity. 

We compute by lattice simulations the derivative F' with respect to the inverse gauge coupling 
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where the subscripts on the averages indicate the value of the external field. Only unconstrained 
plaquette are taken into account in the sum in eq. (2.4). 

In the present work we have considered a static constant abelian chromomagnetic field. The 
continuum gauge potential giving rise to a static constant abelian chromomagnetic field directed 
along spatial direction 3 and direction a in the color space can be written in the following form: 



Af\x) 



(2.5) 



In SU(3) lattice gauge theory the constrained lattice links (see eq. (2.3)) corresponding to the con- 
tinuum gauge potential eq. (2.5) are (choosing a = 3, i.e. abelian chromomagnetic field along 
direction 3 in color space) 
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Since our lattice has the topology of a torus, the magnetic field turns out to be quantized 
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In the following n ext will be used to parameterize the external field strength. 



(2.7) 



3. Deconfinement temperature and critical field strength 

Numerical simulations have been performed on 32 3 x 8 and 64 x 32 2 x 8 lattices, using com- 
puter facilities at the INFN apeNEXT computing center in Rome. A slight modified version of the 
standard HMC R-algorithm [9] has been adopted to simulate QCD with two degenerate flavors of 
staggered fermions of mass am q = 0.075. The critical gauge coupling has been determined from 



the position of the peak in F'[A ext ] (eq. {2A)), the derivative of the free energy with respect to the 
gauge coupling j8, as a function of j8. In figure 1 we show an example of F' measured for « ext = 1 
on a 32 3 x 8 lattice. In the same figure we display also the chiral condensate {yy}. We can see 
that the peak in the derivative of the free energy corresponds to the drop in the chiral condensate, 
the latter being a signal of the transition leading to chiral symmetry restoration. The determination 
of the critical coupling from the peak of the derivative of the free energy is also consistent with 
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Figure 1: The derivative of the free energy eq. ([OJ) with respect to the gauge coupling (left axis, blue 
circles), and the chiral condensate (right axis, red squares) versus /3 . The vertical line represents the position 
of the peak in the derivative of the free energy. 



the determinations obtained by studying the Polyakov and the plaquette susceptibility. Therefore 
we may conclude that the critical coupling of the phase transition can be located by looking at the 
peak of the derivative of the free energy and, as in the case of zero external field and within sta- 
tistical uncertainties, a single transition seems to be present where both deconfmement and chiral 
symmetry restoration take place. 

We have then varied the strength of the external field by tuning up the parameter n ext and 
searched for the phase transition signalled by the peak of the derivative of the free energy finding 
that the critical coupling shifts towards lower values by increasing the external field strength. On the 
other hand we have found that an abelian monopole background field does not have any influence 
on the critical coupling. Indeed we have performed numerical simulations in presence of an abelian 
monopole background field with monopole charge n mon = 10 (again for 2 staggered flavors QCD of 
mass am q = 0.075). The critical coupling has been estimated to be j8 f = 5.4873(192) consistently 
with the critical coupling j8 f = 5.495(25) obtained without any external field. 

Once the critical couplings in correspondence of several values of the external abelian chro- 
momagnetic field have been obtained the corresponding critical temperature is given by 

T c = (R l w , (3.1) 

where L t is the lattice temporal size and a(p c ,m q ) is the lattice spacing at the given critical coupling 
j8 c . We have to face now the problem of fixing the physical scale. In order to reduce the systematic 
effects involved in this procedure we consider the ratios 

T c (gH) _ 
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Figure 2: The critical temperature T c (gH) at a given strength of the chromomagnetic background field 
in units of the critical temperature T c without external field versus the square root of the strength of the 
background field in the same units. Red circles are obtained by adopting the improved scaling function. The 
blue line is the linear best fit. The blue circle on the horizontal axis is the linear extrapolated values for 
the critical background field. Green squares are obtained by adopting the 2-loop scaling function. The blue 
dashed line is the linear best fit. The blue square on the horizontal axis is the linear extrapolated values for 
the critical field. 



where T c is the critical temperature without external field. The above quantities can be obtained 
once the ratio of the lattice spacings at the respective couplings is known. A rough estimate of this 
ratio can be inferred by using the 2-loop scaling function f(g 2 ) 



f(g 2 ) = (bog 2 )-" 1 ^ exp 
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N c is the number of colors and Nf is the number of flavors. A better estimate could be obtained by 
exploiting an improved scaling function f(g 2 )(l +C2a(g) 2 + c^d{g) A ). We do not know, however, 
the values of C2 and c\ for Nf = 2. In a first approximation we will fix C2 and c\ to their quenched 
values [10,11]. 

In figure || the quantities reported in eq. (3.2) are displayed for both choices described above, 
i.e. 2-loop asymptotic scaling and improved scaling. Our main result is that, even in presence of 
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Figure 3: The chiral condensate versus j3 in correspondence of some values of the constant chromomagnetic 
background field. In the inset the region corresponding to the phase transition has been magnified. 



dynamical quarks, the critical temperature decreases with the strength of the chromomagnetic field. 
Furthermore a linear fit to our data can be extrapolated to very low temperatures leading to the pre- 
diction of a critical field strength above which strongly interacting matter should be deconfined at 
all temperatures. Nevertheless, as can be appreciated from figure ||, the exact value of the critical 
field strength is largely dependent on the choice of the physical scale: assuming the 2-loop scal- 
ing function we obtain ^gH c /T c = 26.8(5), while using the improved scaling function we obtain 
yJgH c /T c == 4.29(10); if the deconfinement temperature at zero field strength is taken to be of the 
order of 170 MeV, that means ^gH c in the range 0.7 — 4.5 GeV. Hence in order to get a reliable 
estimate of the critical field strength requires to get a more reliable estimate of the physical scale 
of our lattices. 

To conclude this section we consider our measurements of the chiral condensate. In figure || 
we display the chiral condensate versus the gauge coupling in correspondence of some values of 
the external field strength. Numerical data show that, at least in the critical region, the value of the 
chiral condensate depends on the strength of the applied field. Similar results have been found in 
ref. [12]. 



4. Summary and Conclusions 

We have studied how a constant chromomagnetic field influences the QCD dynamics. By fo- 
cusing on the finite temperature theory we have found that, analogously to what happens in the 
pure gauge theory [2], the critical temperature depends on the strength of the constant chromomag- 
netic background field. More specifically, the critical temperature decreases as the external field 
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is increased and, as an extrapolation of our results, for strong enough field strengths the system is 
always deconfined. A rough estimate of this critical field strength turns out to be of the order of 
1 GeV, which is a typical QCD scale [13]. 

What is more, by comparing the critical couplings determined from the derivative of the free 
energy functional with those determined from the susceptibility of the chiral condensate and of the 
Polyakov loop we have ascertained that, even in presence of an external chromomagnetic back- 
ground field and at least up to the field strengths explored in the present work, the critical temper- 
atures where deconfinement and chiral symmetry restoration take place coincide within errors. 

Another intriguing aspect that deserves further studies is the dependence of the chiral conden- 
sate on the chromomagnetic field strength. 
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